#NEXUS

BEGIN NETWORKS;
Network net = ((((B)#H1,E),(C,A)),(#H1,D));
END;
 
BEGIN TREES;
 
Tree geneTree1 = ((C:3,((B:1,D:1):1,A:2):1):1,E:4);
Tree geneTree2 = (B:4,(D:3,(C:2,(A:1,E:1):1):1):1);
Tree geneTree3 = (D:4,(B:3,((C:1,E:1):1,A:2):1):1);
Tree geneTree4 = (D:3,((B:1,E:1):1,(C:1,A:1):1):1);
 
END;
 
 
BEGIN PHYLONET;
calgtprob net (geneTree1,geneTree2,geneTree3,geneTree4) -bl -o;
END;
===

calgtprob net (geneTree1, geneTree2, geneTree3, geneTree4) -bl -o
Species Network:
((D:1.0,(B:0.3333333333333333)#H1:0.6666666666666667::0.6761569597984691):1.0,((E:0.9814239699997196,#H1:0.6480906366663863::0.32384304020153093):0.018576030000280408,(A:0.9828573145495751,C:0.9828573145495751):0.017142685450424855):1.0);
Total log probability: -23.06958625687312
===
